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Abstract 

We calculate next-to-next-to-leading order soft-gluon corrections to top quark total 
and differential cross sections in hadron colliders. We increase the accuracy of our previ- 
ous estimates by including additional subleading terms, including next-to-next-to-next-to- 
leading-logarithmic and some virtual terms. We show that the kinematics dependence of 
the cross section vanishes near threshold and is reduced away from it. The factorization 
and renormalization scale dependence of the cross section is also greatly reduced. We 
present results for the top quark total cross sections and transverse momentum distribu- 
tions at the Tevatron and the LHC. 



1 Introduction 



The discovery of the top quark in pp collisions at Run I of the Tevatron in 1995 1J and 
its observation currently at Run II, with expected increases in the accuracy of the top mass 
and cross section measurements, have made theoretical calculations of top production cross 
sections and differential distributions an interesting and topical subject. The latest calculation 
for top hadroproduction includes next-to-next-to-leading-order (NNLO) soft-gluon corrections 
to the double differential cross section [2J E] from threshold resummation techniques jUEJE]. 
Near threshold there is limited phase space for the emission of real gluons so that soft-gluon 
corrections dominate the cross section. 

These soft corrections take the form of logarithms, [ln z (x t h)/^th] + ; with / < In — 1 for the 
order a n s corrections, where x t h is a kinematical variable that measures distance from threshold 
and goes to zero at threshold. NNLO calculations for top quark production have so far been 
done through next-to-next-to-leading-logarithmic (NNLL) accuracy, i.e. including leading log- 
arithms (LL) with / = 3, next-to-leading logarithms (NLL) with I = 2, and NNLL with 1 = 1, 
at NNLO Oinj. This NNLO-NNLL calculation has had great success in significantly reducing 
the factorization/renormalization scale dependence of the cross section. Indeed the scale de- 
pendence of top production is almost negligible. However, the dependence of the corrections 
on the kinematics choice is substantial. In Ref. [3J, the top cross section was studied in both 
single-particle-inclusive (1PI) and pair-invariant-mass (PIM) kinematics. Important differences 
between the two kinematics choices were found in both the parton-level and hadron-level cross 
sections, even near threshold. Similar kinematics effects were found for bottom and charm 
hadroproduction j3J[7]. Thus subleading, beyond NNLL, contributions can still have an impact 
on the cross section. If all the NNLO soft corrections are included, there should be no difference 
between the two kinematics near threshold. If all NNLO corrections, both soft and hard, were 
known, there should be no difference between the two kinematics, even far from threshold. 
Away from threshold, where the approximations of Ref. [3] are not expected to apply since 
real emission of hard gluons comes into play, the discrepancy between the 1PI and PIM results 
is not surprising. However, the NNLO-NNLL calculation exhibits some notable discrepancies 
between the two kinematics even at the lowest rj, where rj = s/(4m 2 ) — 1 — > at threshold. 
Thus, additional subleading terms are clearly needed to bring the calculation under further 
theoretical control. 

In this paper, we include additional subleading NNLO soft corrections, including next- 
to-next-to-next-to-leading logarithms (NNNLL) with I = 0, as well as some virtual 5(xth) 
corrections. We apply the method and results of Ref. jH] , based on earlier resummation studies 
jUEHSj, where master formulas are given for the NNLO soft and virtual corrections for processes 
in hadron-hadron and lepton-hadron collisions. As we will see, the subleading corrections do 
indeed bring the 1PI and PIM results into agreement near threshold for both the qq — > tt 
and the gg — ► tt channels, while the discrepancies away from threshold are also diminished, 
especially in the gg channel. Thus the threshold region is brought under theoretical control. 

Since the resummation formalism has been reviewed extensively in Refs. j2 El HI HI IE] , we 
only provide a rough outline here. Threshold resummation is a method of formally calculating 
contributions from soft-gluon emission to all orders in perturbation theory. The resummation is 
normally carried out in moment space where N is the variable conjugate to x t h and the leading 
threshold logarithms are of the form ln 2n N for the order a™ corrections. The resummed cross 
section, in moment space, can then be expanded to NNLO (and even higher orders [2]) and 
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the finite-order result finally inverted back to momentum space. The previous calculations of 
Refs. (21 El and the universal results of Ref. [S] employ this approach. 

In the following section we give the analytical form of the NNLO soft and (some) virtual 
corrections in the qq — ► it channel in both 1PI and PIM kinematics. In Section 3 we give 
the results for the gg — > it channel. Note that while we refer only to it here, the results in 
Sections 2 and 3 are equally valid for all heavy quarks. Section 4 discusses the partonic cross 
sections in both channels. In Section 5 we present the hadronic cross sections and transverse 
momentum distributions for top production in Tevatron Run I and Run II as well as at the 
LHC We conclude with a summary in Section 6. 



2 NNLO soft corrections to qq — > tt 

We first study the qq — > it process. Before discussing the corrections, we introduce our kine- 
matics notation (the same notation will also be used for gg — > it). In 1PI kinematics the process 
is 

q(Pa)+q(Pb) — t(pi)+Xp](p2) (2.1) 

where t is an identified top quark of mass m and X\t\ is the remaining final state that contains 
the?. We define the kinematical invariants s = (p a +Pb) 2 , h = (pb—pi) 2 —m 2 , U\ = (p a —Pi) 2 — m2 
and S4 = s + ti +u\. At threshold, S4 — > 0, and the soft corrections appear as [ln z (s 4 /m 2 )/s4] + . 
In PIM kinematics, we have instead 

q(Pa)+q(p b ) —*it(p)+X(k). (2.2) 

At partonic threshold, s = M 2 , M 2 is the pair mass squared, and t\ = — (M 2 /2)(l — (3m cos#), 

U\ = — (M 2 /2)(l + (3 M cos9) where (3m = \Jl — 4m 2 /M 2 and 9 is the scattering angle in 

the parton center-of-mass frame. The soft corrections appear as [W(l - z)/(l - z)} + with 
z = M 2 /s — > 1 at threshold. 

A more detailed discussion of the kinematics can be found in Ref. 



2.1 qq — > tt channel in 1PI kinematics 

We begin our study with the next-to-leading order (NLO) corrections. In the MS scheme, the 
NLO soft and virtual corrections for qq —>■ it in 1PI kinematics can be written as 
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Here, the Born term is 
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where /j,r is the renormalization scale, Cp = (N 2 — 1)/(2N C ) with N c = 3 the number of colors, 
and K q q = N~ 2 is a color average factor. 
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We also have = ACp and 
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where pip is the factorization scale, Ca = N c , L'^ = [(1 — 2m J s) / 0\ ln[(l — + /?)] and 

/3 = ^/l — Am 2 / s. For later use, we write 



,ipi 



iipi 



2 qq 2 qq 



2C F ln^ N 



so that is the scale-independent part of Cj£L Finally, 



,ipi 



99' 

(1)S+V 1PI 
qqS 



a 



(2.6) 



(2.7) 



where cr^f +V 1PI denotes the <5(s 4 ) terms in Eq. (4.7) of Ref. [Oj with the definitions of t\ and 
Mi interchanged with respect to that reference. We also write 
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where has no scale dependence. 

Before presenting the NNLO soft corrections, we define the constants £2 = 7r 2 /6, £ 4 = 7r 4 /90, 
and (3 = 1.2020569 • • •, and the two-loop constant K = Ca(67/18 — 7r 2 /6) — 5n//9, with rif the 
number of light quark flavors. Finally, we define 

fi{a.) = fi ^ = + (2.9) 



where (3q = (11C^ — 2n/)/3 and 
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Following Ref. [H] we write the NNLO soft-plus-virtual corrections in 1PI kinematics as 
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Here 
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denotes a set of two-loop contributions that are universal for processes with qq initial states 
|S]. Process-dependent two-loop corrections JTD] are not included in but, as we will see 
in Section 4, their contribution is expected to be negligible. The virtual contribution -R^ 1 
is not fully known. However, we can determine certain terms in -R^ 1 exactly. These exact 
terms involve the factorization and renormalization scales as well as the those terms that arise 
from the inversion from moment to momentum space (for a detailed discussion of the inversion 
procedure see Section IIIC and Appendix A of Ref. j2]). 

The terms multiplying ^(s^) involving the factorization and renormalization scales are given 
explicitly by 



7 B IPI 

qq 



a 2 (/i 2 



7P 



171 



c F 

2 



In 



tlU! 



rrr 



2 



rrr 



+ In 
+ In 



m 2 
m 2 



/4 



hi 



In 

-2C F ( 2 
C F 



4 



In 



ti Ui 



rrr 



rrr 



f In' 



3/3 2 



16 



In 



rrr 



T. 



ipi 

2qq 



2C F ln 



m 
s 



K 



SC^Ca + CVyln 



3/?o 
4 



In 



m 4 



m 



m 



-o H— +^ln — +Ti 



rrr 

ipi 

99 



where 



/(2) 



( 32 - I C2 + ^ 3 



11 



17 



12 



96 



C2 



7C3 + — C2 + — + n f C F -- - — 



6 48 



(2.13) 



(2.14) 



The terms multiplying 5(s^) resulting from inversion (( terms) that do not involve the 
factorization and renormalization scales are (21 IE] 
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2.2 



qq 



tt channel in PIM kinematics 



Next, we study the soft-gluon corrections in PIM kinematics. The MS NLO soft and virtual 
corrections to qq — > tt in PIM kinematics are 
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where |pim indicates that for ti, U\ we use the expressions below Eq. ()2.2j) . Also c PI1 ^ = 4CV 
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and 
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Note that the scale-independent Tfg- 1 is related to its 1PI counterpart by 
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Here a ^ S+MF denotes the soft and mass factorization subtraction terms calculated in Ref . [3] . 
The prime indicates that we drop the overall 5(1 — z) or £(34) coefficient from the expressions 
in Eqs. (82), (A8), and (A9) of Ref. [3]. 

In PIM kinematics, the NNLO soft-plus-virtual corrections are 
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Again, only certain terms in R^™ that can be determined exactly are included. The terms 



multiplying 5(1 — z) that involve the factorization and renormalization scales are 
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The terms multiplying 8(1 — z) that arise from inversion and do not involve the factorization 
and renormalization scales are given by 
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3 NNLO soft corrections to gg 
3.1 gg — > channel in 1PI kinematics 



We now turn to the gg channel. We write the MS NLO soft-plus-virtual corrections for gg — > tt 
in 1PI kinematics as 
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The Born term is given by 
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Finally, 
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where here cr^Js +V1FI denotes the scale-independent <5(s 4 ) terms in the NLO cross section. 
These terms are given by Eq. (6.19) in Ref. 

The NNLO soft-plus-virtual corrections in 1PI kinematics are 



2 d 2 a$ 1PI 
-i gg _ p 



dti du 



yy 



71^ 



o ( C 3 gg) 



+ 



+ 



IPI 1PI , ( r lPl\ 2 /- ( r 
c 3gg c lgg + \ c 2gg) ~ ^2 [C, 



1PI 

3 gg 



ln 3 (s 4 /m 2 ) 

s 4 

/^Orpipi , A 



+ 

— i — c lpI 

2 2 99 1 4 3 gg 



3 ipi ipi _ A) ipi 

2 C 3 gg C 2 ^ c 3 gg 



r ln 2 (s 4 /m 2 ) 



s 4 







ln(s 4 /m 2 ) 






s 4 



ipi ipi _ > ipi ipi , > opi \ 2 , Ao ipi , 

C 2gg C lgg S,2 C 2 C3 gg + C, 3 ^Cg gg J -T C^ gg m 



Po 1 2 ( V^F 



/4 



+ CU^fln 2 — -C A AT In ^ — CaKIu 



m 4 



+ 0, 



(2) 



+ 



4f ,,2 \ 



a 4 (/i 



-r lpI 4 C 
2 3 gg ^gg 



ln 2 (s 4 /m 2 ) 



2 cjw -^U 



A) 



IPI rpdPI 



^2gg 



.'/.'/ 



99^ J- 99 1 99 



C 3 gg J 1 



ln(s 4 /m 2 
s 4 



+ 



IPI J- IPI , r-u 1 
C l gg - C2C3 M H~ — In 



Ao 
4 



^gg + I c 



.in _ A) 

2gg 2 







■ 1 ■ 






-S4- 



where 



(3.8) 



F c 
gg 



71 



-K gg B QED (N'-l) 2 In -i 



wA (t 2 - u 2 ) 



ti 



4 rff + 2(A^ 2 -2)rff 



1 - 
4 



1t\U\ 



4(rff) 2 + (iV 2 + 4) In 2 (| 



(rff) 2 - 2(rfi) 5 



2i\Un 2 (Hi 
£1 



(3.9) 



with 



rff = -c F (L'p + i)+c A 



^yy 
22 



2 + ln 



Here 



(3.10) 

, ( , . , (3.11) 

2 3 ' 108/ 1 ° r 1 ' 

denotes a set of universal two-loop contributions for processes with gg initial states process- 
dependent two-loop corrections [TO] are not included in Qfg- 



8 



The contributions R^ 1 and -R^ PI are virtual corrections that are not fully known. As in the 
qq channel, we keep only certain terms that can be determined exactly. The terms multiplying 
5(s4) involving the factorization and renormalization scales are 



7 B IPI 

99 



L2 (»f\ \ C a ,„ 2 5f3 °-C A \n ( tlUl ) + 3/9 o 



In 



W 



rrr 



rrr 



rrr 



16 



2C2C. 



+ ln(4V^- 



G A In 



*lUi 

m 4 



2 



In 



2 ( f J, R\ 3/^o 



16 



C\ In 2 ( ^ In ( ^ ) - In f ^ In f ^ 



+ 2C A ( 2 



rrr 



T? g ]-2C A \n 



rrr 



m 



rrr 



+ 



IT 2 



2C A ( 2 A c gg + (c A ln( 



h Mi 

m 4 



7? 



clPI 

93 



\m 2 4 



3/5 °T clpI In 

1 l99 m 



'/4 N 

m 2 



where 



/(2) ^2 / 2 , 3 a \ (C F C A \ 



(3.12) 
(3.13) 



The terms multiplying 5{sa) that arise from inversion and do not involve the factorization 
and renormalization scales are 
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3.2 — > tt channel in PIM kinematics 

We continue our study of subleading terms in the gg channel by writing the MS NLO soft-plus- 
virtual corrections for gg — > it in PIM kinematics as 
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Here the Born term is 
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Here a' S+MF denotes the soft and mass factorization subtraction terms calculated in Ref . |3j ■ 
The prime indicates that we drop the overall 5(1 — z) or 5(s$) coefficients from the expressions 
in Eqs. (82), (A10), and (All) of Ref. 0. 

The NNLO soft-plus-virtual corrections in PIM kinematics are 
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where A c gg and F gg are the 1PI functions given in the previous subsection. 

The virtual corrections -Rg<j M , Rgg are a l so n °t fully known. We keep only certain terms 
that are determined exactly. The terms multiplying 5(1 — z) that involve the factorization and 
renormalization scales are 
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The terms multiplying 5(1 — z) that arise from inversion and do not involve the factorization 
and renormalization scales are 
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4 Partonic cross sections 

Any difference in the integrated cross sections due to kinematics choice arises from uncalculated 
subleading terms. At leading order (LO) the partonic threshold condition is exact and there is 
no difference between the total cross sections in the two kinematic schemes. However, beyond 
LO additional soft partons are produced and there is a difference when not all terms are known. 
The total partonic cross section may be expressed in terms of dimensionless scaling functions 
flj that depend only on rj = s/Am 2 — 1 

t?) = (W/i^E $%) In' (4) • (4.1) 

m k=o 1=0 \ m J 

These scaling functions all multiply powers of ln(/i 2 /m 2 ) and thus do not depend on \x them- 
selves. Here we have set /x = \ip = [Mr. We work in the MS scheme throughout. 

Previously, we constructed LL, NLL, and NNLL approximations to fjj in the qq and gg 

channels for k < 2, I < k jSJ- We now present the full soft-plus- virtual results for the fjf and 

(2 2) (2 0) 

' scaling functions and the partial results for ' that include the soft NNNLL and those 
virtual terms calculated in sections 2 and 3. 

(2 1) 

We begin with a comparison of the full soft-plus-virtual 1PI and PIM contributions to ' 

and fif' 2 \ shown in Fig. [TJ The upper plots are for the qq channel. The left-hand side of Fig. [T] 

(2 1) 

compares the 1PI and PIM scaling functions for f^' . At low 77, closer to partonic threshold, 

the agreement is very good, better than that obtained at NNLL in Ref. 0. The agreement 

(2 2) 

is also improved at large r/. The right-hand side shows the / ? y scaling functions in both 

kinematics. The results for f^ 2 ^ 1 remain unchanged from those of Ref. [Sj. The agreement 
between the two kinematics choices is excellent. 

The lower plots of Fig. ^ show the corresponding scaling functions in the gg channel. The 
agreement between the two kinematics choices is somewhat improved at high 77 as compared to 
previous NNLL results [3J . We note that there is some ambiguity in the way that the expressions 
for the gg partonic cross sections can be written at threshold. We have investigated the effect 
of replacing 1 — 1t\U\j s 1 with {t\ + u^)/s 2 in Eq. ()3.6j) . more consistent with the expressions in 
Ref. [TT]. These two expressions are equivalent at threshold, S4 = and z — 1, but can differ 
at large rj. Note that fg 2 ,' 2 ^ is not affected by this replacement. The resulting differences in 
fg 2 ,' 1 ^ are small, appearing only at i] > 0.1 where the agreement between the scaling functions 
in the two kinematics begins to diverge. The main effect of the second choice is to make the 
PIM result for fgg' 1 ^ more negative at large rj. We thus use the expressions as written in the 
text to be consistent with those of Ref. jS] . 

We now turn to the fj; 2 ' ^ scaling functions, the most important contributions at NNLO and 
independent of ln(/i 2 /m 2 ). We add the NNNLL terms, i.e. terms proportional to [1/34]+ (1PI) 
and [1/(1 — z)] + (PIM), to our previous NNLO-NNLL calculation. We also investigate the 
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Figure 1: The MS scheme scaling functions multiplying the scale-dependent logarithms, 

(left-hand side) and f^' (right-hand side). The upper plots are for the qq channel while the 
lower plots are for the gg channel. The solid curves are for 1PI kinematics, the dashed for PIM 
kinematics. 



effect of keeping the virtual ( terms resulting from the inversion from moment to momentum 
space. To demonstrate the effect of adding successive subleading contributions, in Fig. |21 we 
show the NNLL results in the upper plots, the scaling functions through NNNLL in the middle 
plots, and the results with the NNNLL and virtual ( terms in the lower plots. 

We first discuss the results for the qq channel in the MS scheme, shown on the left-hand 
side of Fig. El Note that to NNLL, the two kinematics choices give rather different results, even 
at low r\. When the NNNLL contribution is added, both the 1PI and PIM results are reduced 
relative to the NNLL over all rj. The agreement between the two kinematics is much improved 
up to rj > 0.01. Adding the virtual ( terms resulting from inversion improves the agreement 
between the 1PI and PIM kinematics further for 0.01 < i] < 0.1. At rj > 0.1, the region where 
the parton luminosity peaks for tt production at the Tevatron, the additional virtual ( terms 
provide a further small reduction. With the subleading terms, the 1PI result is smaller than 
previously but positive while the PIM result becomes more negative. However, on the whole, 
the subleading terms bring the 1PI and PIM results into better agreement over all r]. We note 
here that the effect of the virtual ( terms is numerically small, as is also the case for the gg 
channel and for the hadronic results for both channels in the next section. This small effect is in 
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Figure 2: The f-j ' scaling functions in the MS scheme. The left-hand side shows the results 
for the qq channel while the right-hand side shows the results for the gg channel. The top plots 
show the NNLL result from Ref. The center plots give the results through NNNLL and 
the bottom plots give the results including the virtual ( terms. The solid curves are for 1PI 
kinematics, the dashed for PIM kinematics. 



agreement with the arguments in Section IIIC of Ref. [2] concerning resummation prescriptions. 
There it was shown that when subleading terms from inversion are calculated exactly they do 
not have an unwarrantedly large effect on the numerical results. 

A similar trend is seen for the gg channel on the right-hand side of Fig. The agreement 
between the NNLL 1PI and PIM scaling functions at low 77 is significantly better than in the 
qq channel. This may perhaps be a consequence of the more complex color structure of the gg 
channel. Note however the significant divergence at large 77. The 1PI NNLL result is large and 
positive while the PIM is large and negative. Again, inclusion of the subleading contributions 
improves agreement over all rj. There is only a small improvement possible at low 77. However, 
the improvement at larger rj, rj > 0.1 is notable. The 1PI result with soft NNNLL plus virtual 
( terms is reduced by nearly a factor of two relative to the NNLL result at 77 = 10. Likewise, 
the subleading terms stop and reverse the downward trend of the PIM scaling functions. The 
1PI gg contribution will still be positive while the PIM will still be negative but the difference 
may not be as large as before. Using the alternate expression, {t\ + uD/s 2 , in Eq. flUS)) does 
not significantly change the results, particularly for 1PI kinematics. The PIM result becomes 
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slightly more negative at intermediate rj, r\ « 1. 

Finally we note that if we had kept only the ( contributions in the [l/s 4 ] + and [1/(1 — z)] + 
terms the 1PI and PIM results would not have agreed near threshold. The full NNNLL result, 
given in sections 2 and 3, is required for the result to be independent of kinematics choice near 
threshold. This agreement also indicates that additional two-loop contributions not included 
in our expressions should be small. 

We now turn to our calculations of the hadronic total cross sections and transverse momen- 
tum distributions. 



5 Hadronic total cross sections and £>r distributions 

The inclusive hadronic cross section is obtained by convoluting the inclusive partonic cross 
sections with the parton luminosity, defined as 



L 1 

(r,/4) = T J dx 1 Jdx 2 5(xiX2-T)(j) i / hl ( y x 1 ,fi 2 F )(l) j / h2 (x2,iJ 2 F ), (5.1) 



where (pi/h{x, fip) is the density of partons of flavor i in hadron h carrying a fraction x of the 
initial hadron momentum, at factorization scale /ip. Then 



dr 



/(IT 
— $ y (r,/4) (Tij(TS,m 2 ,fi 2 F ) (5.2) 

*,3=q,q,g 4m 2 /s T 

E / d\og 10 r]—— ln(10)$ij(77, fj.%) 0^(77, m 2 , /4) 

.. .. - _ J —00 1 T" <] 



where 



^ 4m 2 4m 2 ' ^ 

and S is the hadronic Mandelstam invariant. Our investigations in Ref. [3] showed that the 
approximation should hold if the convolution of the parton densities is not very sensitive to the 
high rj region. 

We use the recent MRST2002 NNLO (approximate) parton densities [12 with an NNLO 
evaluation of a s . The parton luminosities, weighted to emphasize the most important contribu- 
tions to the hadronic cross sections, are shown for v^S* = 1.96 TeV in Fig. The qq luminosity 
is nearly 50% higher than the CTEQ5M [inj qq luminosity used in Ref. [Sj. The gg luminosities 
for the two sets are rather similar. The peak of the luminosity is at rj < 1, but still in a regime 
where the 1PI and PIM results differ most. Fortunately the gg luminosity is small compared 
to the qq luminosity since the differences in the kinematics is largest in the gg channel. 

Our calculations use the exact LO and NLO cross sections with the soft NNNLL and virtual 
£ corrections and the full soft-plus- virtual scale-dependent terms at NNLO. In addition we 
multiply the NNLO scaling functions by a damping factor, + rj, as in Ref. to lessen 

the influence of the large r] region where the threshold approximation does not hold so well. 

In Fig. 0J we present the NLO and approximate NNLO it cross sections at y/S = 1.8 TeV 
(left-hand side) and 1.96 TeV (right-hand side) as functions of top quark mass for \i = m. The 
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Figure 3: The qq (left-hand side) and gg (right-hand side) parton luminosities in pp collisions 
at = 1.96 TeV. The solid curves are calculated with ji — m — 175 GeV while the upper 
dashed curves are with /1 = m/2 and the lower dashed curves with /1 = 2m. 



NNLO results include the soft NNNLL and virtual ( terms in 1PI and PIM kinematics. We 
also show the average of the two kinematics results which may perhaps be closer to the full 
NNLO result. Here the NNLO PIM cross section is slightly lower than the NLO cross section 
for all masses shown. In Ref. the PIM cross section was a bit higher than the NLO. The 
reduction of the PIM qq result, dominant for pp — > it, lowers the total PIM cross section. The 
NNLO 1PI cross section remains above the NLO for all m although the NNLO cross section is 
not as large as previously, due to the subleading terms. The average of the two kinematics is 
just above the NLO cross sections for both energies. 

Going to higher scales increases all the NNLO corrections so that both kinematics choices 
give cross sections larger than the NLO. On the other hand, at lower scales, the NNLO cross 
sections are reduced relative to the NLO. The ratio of the NNLO to the NLO cross sections, 
the K factors, are shown in Fig. 03 as functions of mass for fi = m (upper plot), 2m (middle 
plot and m/2 (lower plot) at v^S* = 1.8 TeV. In keeping with the results in Fig. 01 for fi = m 
K < 1 for PIM kinematics, > 1 for 1PI and for the average. The K factors are larger for 
fi = 2m and smaller for fi = m/2. Note also that K is almost independent of m. We remark 
that the NLO/LO K factor, while also essentially mass independent, is typically larger than 
the NNLO/NLO K factors shown here. It is ~ 1.25 for fi = m, 1.52 for fi = 2m and 0.94 
for n = m/2. Only the last value is similar to that of the NNLO/NLO average K factor in 
Fig. E| The small K factors, obtained with results calculated with the MRST NNLO parton 
distribution functions at each order, indicate good convergence. Even though the results are 
shown at ^/S = 1.8 TeV, the K factors at = 1.96 TeV are very similar. 

We now examine the scale dependence in Fig. as a function of top quark mass and in 
Fig. El as a function of fi/m with m = 175 GeV. Figure shows the ratio of the cross sections 
with fi = 2m to fi = m on the left-hand side and the ratio for fi = m/2 to fi = m on the 
right-hand side at both NLO and NNLO at v^S* = 1.8 TeV. The ratios are nearly independent 
of mass at this energy. The scale dependence is reduced at NNLO relative to NLO. The NNLO 
results are very similar for the two ratios. In contrast, the LO scale dependence is much larger, 
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Figure 4: The it total cross sections in pp collisions at y S = 1.8 TeV (left-hand side) and 1.96 
TeV (right-hand side) as functions of m for \i — m. The NLO (solid), and approximate NNLO 
1PI (dashed), PIM (dot-dashed) and average (dotted) results are shown. 



cr(/i = 2m)/ff(/i = m) « 0.74 and a(fi = m/2)/cx(/z = m) « 1.4. The difference between the 
scale dependence at \/S =1.8 TeV and 1.96 TeV is negligible. 

We have also calculated the cross sections as functions of fi/m for 0.2 < \ijm < 10 at 
\/S = 1.96 TeV and m = 175 GeV in Fig. [7| The NLO cross section is not as strong a 
function of fi/m as the LO cross section. In fact, it is seen to rise with n/m and then turn 
over. The NNLO cross sections, however, exhibit even less dependence on fi/m, approaching 
the independence of scale corresponding to a true physical cross section. They change by less 
than 15% over the entire range of \ijm considered. The change in the NNLO cross sections 
through the range m/2 < /i < 2m, normally displayed as a measure of uncertainty from scale 
variation, is less than 3%. Note also that, at this energy, the absolute difference between the 
1PI and PIM cross sections is also not large. 

In Table 1, we give the NLO, NNLO 1PI and NNLO PIM it total cross sections at = 1.8 
and 1.96 TeV for pp interactions, corresponding to Tevatron Runs I and II. The results are 
presented for m = 175 GeV and fi = m/2, m, and 2m. We show the results of our calculations 
with the MRST2002 NNLO parton densities ^2] and the three-loop a s . We compare these 
with results of calculations with the CTEQ6M NLO parton densities [13] and the two-loop a s . 
The results with the two different sets of parton densities are quite similar even though the 
parton densities are evaluated to different orders. Note that the NNLO scale dependence is 
negligible compared to the NLO scale dependence. The kinematics dependence of the NNLO 
cross sections thus remains the largest source of uncertainty. At y/S = 1.8 TeV, averaging over 
the 1PI and PIM NNLO results with the two sets of parton distributions at /i = m = 175 GeV, 
our best estimate for the cross section is 5.24 ± 0.31 pb where the quoted uncertainty is from 
the kinematics dependence. At y/S =1.96 TeV our corresponding best estimate is 6.77 ± 0.42 
pb. 

We note that the cross sections presented in Table 1 are significantly lower than our previous 
estimates [21 IS] at both NLO and NNLO. The difference at NLO is solely due to the new sets of 
parton densities used here, MRST2002 and CTEQ6M, relative to CTEQ5M in Refs. HE]. Our 
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Figure 5: The NNLO K factors at y S = 1.8 TeV as functions of top quark mass in pp collisions 
with fi = m (upper), \x = 2m (middle) and // = m/2 (lower). The curves show the ratio of the 
approximate NNLO 1PI (dashed), PIM (dot-dashed) and average (solid) cross sections to the 
NLO cross section. 



new NLO results are around 3% lower. The effect of the new densities on the NNLO corrections 
is even larger. The NNLO-NNLL 1PI corrections are smaller than our previous results [21 El 
by 14% for fi = m and 18% for fi = 2m with the MRST2002 NNLO densities. Most of this 
difference is in the relative values of a s between the two sets. In addition, the new subleading 
terms we have included here further reduce the magnitude of the NNLO corrections. The 
combined effect of the new parton densities and new subleading terms make our new estimates 
for the total NNLO ti cross section noticeably smaller. 

In Fig. |H1 we show the top quark transverse momentum distributions at y~S = 1.8 and 
1.96 TeV. The NLO and NNLO 1PI results are shown using the MRST2002 NNLO densities. 
Details of the hadronic calculation of the px dependence are given in Appendix B of Ref. j2|. 
At NNLO we observe an enhancement of the NLO distribution with no significant change in 
shape. This pattern agrees with earlier resummed results on top transverse momentum and 
rapidity distributions [21 HH] ■ 

Finally we discuss top production in pp collisions at the LHC The weighted parton lumi- 
nosities are shown in Fig.[U]for the maximum LHC pp energy, v^S* = 14 TeV. The gg luminosity 
now dominates the qq by a factor of 4. The peak of the luminosity is still at rj < 1 so that this 
energy is not very far from partonic threshold. However, large uncertainties may be expected 
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Figure 6: The scale dependence of the it total cross sections in pp collisions at v^S* = 1.8 TeV 
as a function of top quark mass. The left-hand side shows the ratio (fi = 2m) /(/i = m) while 
the right-hand side gives the ratio for (fi = m/2)/(p = m). The NLO (solid), and approximate 
NNLO 1PI (dashed), PIM (dot-dashed) and average (dotted) results are shown. 



in the gg channel since the difference in the kinematics choice, largest in this channel, will be 
emphasized by the high gg luminosity. 

Since the gg contribution dominates at high energy, the difference in the results for the two 
kinematics increases strongly with energy. The complex color structure of the gg channel may 
be better suited to 1PI kinematics and thus this kinematics choice could be more appropriate in 
processes where the gg channel dominates, see Ref. |7] for discussions of bottom and charm pro- 
duction. The NNLO 1PI scale dependence at high energy seems to support such a conclusion. 
At v^S" = 14 TeV, the NNLO 1PI scale dependence is 4%, smaller than the 9% dependence of 
the NLO cross section, an acceptable behavior, similar to that at the Tevatron. However, the 
NNLO gg PIM contribution is large and negative. The qq PIM contribution is also negative for 
fi < m albeit much smaller than the gg contribution. The NNLO PIM cross section is reduced 
by nearly a factor of two relative to the NLO. The scale dependence is similarly large. Thus 
we will only provide NNLO 1PI results for the LHC. At y/S = 14 TeV with m = 175 GeV and 
the MRST2002 NNLO parton densities, the NLO cross section is 808.8 pb for y. = m/2, 794.1 
pb for \i — m, and 744.4 pb for \i = 2m. The corresponding NNLO 1PI cross sections are 845.2 
pb for jj, = m/2, 872.8 pb for \x = m, and 875.1 pb for y = 2m. In Fig. Elwe show the NLO 
and NNLO 1PI top quark pr distributions at v^S* = 14 TeV. Here also the NNLO corrections 
enhance the NLO result without a change in shape. 



6 Conclusions 

In this paper we have calculated soft NNLO corrections to the total top quark cross section and 
top transverse momentum distributions in hadron-hadron collisions. We have added new soft 
NNNLL terms and some virtual terms, including all soft-plus-virtual factorization and renor- 
malization scale terms. We have found that these new subleading corrections greatly diminish 
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Figure 7: The scale dependence of the it total cross sections in pp collisions at VS = 1.96 TeV 
as a function of /i/m. The LO (dot-dot-dot-dashed), NLO (solid), and approximate NNLO 1PI 
(dashed), PIM (dot-dashed) and average (dotted) results are shown. 

the dependence of the cross section on the kinematics and on the factorization/renormalization 
scales. We have provided numerical results for the total cross section and top transverse mo- 
mentum distributions for top quark production at the Tevatron, at both Run I and II, and at 
the LHC. 
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Figure 8: The top quark transverse momentum distribution with m = 175 GeV at v^S* = 1.8 
TeV (left) and 1.96 TeV (right). The NLO (solid /i — m, dotted \i = m/2, dot-dashed \x = 2m), 
and approximate NNLO (/i = m) 1PI (dashed) results are shown. 
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Figure 9: The qq (left-hand side) and gg (right-hand side) parton luminosities in pp collisions 
at V~S = 14 TeV. The solid curves are calculated with // = m — 175 GeV while the upper 
dashed curves are with ^l — m/2 and the lower dashed curves with /j, = 2m. 
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Figure 10: The top quark transverse momentum distribution with to = 175 GeV at = 14 
TeV. The NLO (solid /i — m, dotted /x = m/2, dot-dashed /i = 2m), and approximate NNLO 
(/i = to) 1PI (dashed) results are shown. 
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